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§1. Introduction 

In 1949, a theorem in quantum mechanics was established by Levinson[l], which relates 
the phase shift at threshold (zero momentum) for the Zth partial wave, 5i(0), to the number 
of bound states with the same azimuthal quantum number, m. This is one of the most 
interesting and beautiful results in nonrelativistic quantum theory. The subject was then 
studied by many authors, some are listed in the References [2-5]. The relativistic general- 
ization of Levinson's theorem has also been well established [6-8]. However, most of these 
authors deal with the problem in ordinary three-dimensional space. To our knowledge a 
two-dimensional version of Levinson's theorem was not presented in the literature. The 
purpose of this work is to develop an analogue of this theorem in two spatial dimensions. 
It relates the phase shift at threshold for the rath partial wave, rj m (0), to the total number 
of bound states with angular momentum mh, n m (the total number of bound states with 
angular momentum — mh is also n m ): 

r] m (0)=n m n, 771 = 0,1,2,.... 

This is similar to but slightly simpler than the original one in three dimensions. In three 
dimensions Levinson's theorem takes the same form with m replaced by I. But when 
I = the relation must be modified if there exists a zero-energy resonance (a half-bound 
state). The mathematical origin is that the behaviour of the phase shifts 5i(k) near k — 
may be different for / = and 1^0. In two dimensions no similar situation occurs. 
This paper is organized as follows. In the next section we give a brief formulation of the 
partial-wave method for nonrelativistic scattering in two spatial dimensions. In §3 we 
discuss the behaviour of the phase shifts near k = in some detail. In §4 we establish the 
Levinson theorem using the Green function method[2,3,5]. §5 is devoted to the discussion 
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of some aspects of the theorem. 

§2. Partial-wave method in two dimensions 

A particle with mass ji and energy E moving in an external field V(r) satisfies the sta- 
tionary Schrodinger equation 

f> 2 

Hip = -— V 2 -0 + V(t)i/> = Etj). (1) 

We use the polar coordinates (r, 9) as well as the rectangular coordinates (x, y) in two 
spatial dimensions. For scattering problems E > 0(we assume that V — > more rapidly 
than r~ 2 when r — > oo). The incident wave may be chosen as 

A = e lkx (2) 



which solves Eq.(l) when r — > oo provided that k = \j2/j,E/fi 2 . The scattered wave should 
have the asymptotic form 

^ ^ ^-f{9)^ (3) 
where the factor \f% = e i7F//4 is introduced for later convenience. This also solves (1) when 
r — > oo. It is easy to show that the differential cross section cr(0)(m two spatial dimensions 
the cross section may be more appropriately called cross width) is given in terms of the 
scattering amplitude f(9) by 

(7(0) = |/(0)| 2 . (4) 
The outgoing wave comprising (2) and (3) and thus takes the following form at infinity. 

i> ™ e lhx + \j*-f(6)e lkr . (5) 

We are interested in central or spherically symmetric (actually cylindrically symmetric 
in two dimensions) potentials V(r) = V(r). In this paper we deal only with central 



potentials. Then solutions of Eq.(l) may be expanded as 

+00 

V>M)= E a m R\ m \{r)e me (6) 

m=— 00 

where R m (r) satisfies the radial equation 



R'L+ l -R' m + U 2 -^V-^-)R m = 0, m = 0,l,2,.... (7) 



If V = 0, the regular solution of (7) is the Bessel function and may be taken as R$(r) = 
\fkJ m (kr) and hence has the asymptotic form 

*S?M ™ {l<™ (*r " ™ " j) • (8) 

We have assumed that V(r) — > more rapidly than r~ 2 when r — > 00, so the solution 
R m (r) approaches the form of a linear combination of the Bessel function and the Neumann 
function at large r and may have the asymptotic form 



R m (r) ™ 1/ — cos 



nr 



. mil it 

kr - — - - + Vm{k) 



(9) 



Here r) m (k) is the phase shift of the mth partial wave. It is a function of k. As in three 
dimensions, all r) m (k) are real in a real potential. Substituting (9) into (6) gives one 
asymptotic form for ip, while substitution of the formula 



E i^J\rn\{kry me (10) 



and the asymptotic form of the Bessel functions into (5) gives another. Comparing these 

two asymptotic forms one finds an expression for f{6) in terms of the phase shifts: 

+00 / 9 

W)= E J-re^MsmTjHe*"*. (11) 

m =_oo V 7™ 

The total cross section a t turns out to be 

<7 t = ^ ^ o-(0) = - (^sin 2 1)„ + 2^ sin 2 r) m J . (12) 



From the above relations one easily realizes that all information of the scattering process 
is contained in the phase shifts r) m (k). The latter are determined by solving the radial 
equation (7) with the boundary condition (9) and thus depend on the particular form of 
V(r). In an attractive field, the number of bound states with given angular momentum 
mh, denoted by n m above, also depends on the particular form of V(r). It will be shown 
that n m is related to r) m (k) at threshold. This is similar to Levinson's theorem in three 
dimensions. In the next section we first discuss the behaviour of r) m (k) near k — 0. 

§3. Phase shifts near threshold 

Assuming that V(r) is less singular than r -2 when r — > 0, then the regular solution of the 
radial equation (7) may have the following power dependence on r near r = 0. 

n r m 

fm(r,k) r -^—- m = 0,l,2,.... (13) 
z m m\ 

Here we denote the regular solution of (7) with the boundary condition (13) by f m (r, k). 
Note that the equation (7) depends on k only through k 2 , which is an integral function of 
k, and the boundary condition (13) is independent of k. Thus a theorem of Poincare tells 
us that f m (r, k) is an integral function of k for a fixed r. On the other hand, the solution 
R m (r) with the boundary condition (9), which is proportional to f m (r, k), need not be an 
integral function of k. We denote a potential that satisfies V(r) = when r > a > by 
V a (r). In such potentials the solution of (7) when r > a may take the form 

R+(r) = y/k[cosrj m J m (kr) - sin rj m N m (kr)] : m = 0,1,2,... (14) 

where the superscript "+" indicates r > a. It is easy to verify that indeed satisfies 

the boundary condition (9). When r < a we have 

R-(r)=A m (k)f m (r,k), m = 0,1,2,... (15) 
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where the superscript "— " indicates r < a. In general the coefficient A m depends on k, 
so that the two parts of R m (r) can be connected smoothly at r = a. This leads to 



where p = ka and 



tan, m = P 4¥\ Z ^¥A (16) 
' P N' m {p) - l3 m {p)N m {p) 1 ; 



^> = 7$f (17) 



where the prime indicates differentiation with respect to r. As mentioned above, f m (a, k) 
is an integral function of k, so is f' m {a, k). Moreover, both of them are even functions of 
k since Eq.(7) depends only on k 2 . Therefore when k — > or p — > 0, the leading term for 
Pm(p) ma Y have one of the following forms 

Pm{p) -> a m P~ 2lm , 

where a^, a m , and 7 m are nonzero constants, while 1^ and l m are natural numbers. Using 
these relations and the leading terms of J m (p) and N m (p) for p — > 0, the leading term in 
tanr/ m when p — > can be explicitly worked out. When 7 m = — m some care should be 
taken. However, careful analysis gives in any case 

tan77 m ^fe m p 2 ^ or (k — > 0) (18) 

imp 

where 6 m 7^ is a contant and p m is a natural number. Eq.(18) is important for the 
development of the Levinson theorem in the next section. 

For comparison we give the corresponding results in three dimensions. The phase 
shifts are denoted by Si(k). By similar analysis it can be shown in a potential V a (r) that 

tan<&i->cjp 2?, -\ 1 = 0,1,2,... (k — > 0) (18') 



where c\ ^ is a constant, and q\ is a natural number for Z 7^ 0, while go m ay be a natural 
number or zero. We see that 5/(0) generally equals a multiple of n for all /. But 5 (0) 
gets an additional 7r/2 when q = 0. The latter case does not occur for any rj m (0) : which 
is obvious from (18). The difference between (18) and (18') comes from the fact that the 
Neumann function in the solution (14) involves the logarithmic function while the spherical 
Neumann function in the three-dimensional solution does not. It can be shown that 5 (0) 
gets an additional n/2 when there exists a half-bound state(a zero-energy resonance) in 
the angular momentum channel / = 0. 

§4. The Levinson theorem 

Now we proceed to establish the Levinson theorem by the Green function method. Intro- 
duce the retarded Green function G(r, r', E) defined by 

where {^(r)} is a complete set of orthonormal solutions to (1), and e = + . G(r,r',E) 
satisfies the equation 

(E-H + ie)G(r,r',E) =5(r-r'). (20) 
For a free particle we have a similar definition: 

G<0,(r ' r ' £) =?fSr? (21) 

where {ijj^\r)} is a complete set of orthonormal solutions to Eq.(l) with V = 0. G^r, r', E) 
satisfies 

(E-H + ie)G (0) (r, r', E) = 5(r - r) (22) 
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where H is the Hamiltonian of the free particle. We have the integral equation for 
G(r,v',E): 

G(r,r',E)-G^(r,r',E) = J dr" G (0) (r, r", E)V(t")G(t", r', E). (23) 
In a central field V(r) = K(r)(not necessarily V a (r)), we have 



u\ m \ K (r) e 



imd 



Mr, 0) = ^w(r, 9) = |f "'P ' y= , m = 0, ±1, ±2, . . . 

where n is a quantum number associated with the energy E mK which is determined by 
solving the radial equation 



-£2{E mK -V) — 



u mK = 0, m = 0,l,2,... (24) 



ft r 2 

with appropriate boundary conditions in the radial direction. The radial wave functions 
u mK {j-) satisfy the orthonormal condition 

(umK,u mK >) = dru* mK (r)u mK >(r) = 5 KK >. (25) 
Jo 

For an attractive field we have discrete spectrum(_E mK < 0) as well as continuous spectrum 
(E mK > 0). We can, however, require the wave functions to vanish at a sufficiently large 
radius R(R ^> a for V a (r)) and thus discretize the continuous part of the spectrum. In 
this case the upper limit of the integration in (25) should be replaced by R. It is easy to 
show that 

G(r, r', E) = G (r, r', + £ G m (r, r', (26) 

27r m=l 71 

where 

^^) = Et^T^' m = 0,l,2,.... (27) 
„ Vrr'{E - E mK + ie) 



For a free particle, the following results can be obtained in the same way. 
G<°> (r, r', £) = <#> (r, r', E) ±- + ± Gg) (r, r', £) 

2?I " m=l T 



(28) 



where 

? .(o) ( r ) v (o)*( r >\ 

Gg ) (r,r / ,g)=E Ari y " = 0,1,2,... (29) 

« VrP(E - EinL + ie) 

where u^(r) satisfies (24) with V — 0, and the energy spectrum(£'^ > 0) is discretized 
according to the above described prescription. Thus the orthonormal relation for u$ K (r) 
is similar to (25). Substituting (26) and (28) into (23) we get an integral equation for 
G m (r,r',E): 

G m (r, r', E) - G%\r, r>, E) = J dr" r"G%\r, r", E)V(r")G m (r", r',E), m = 0, 1, 2, . . . . 

(30) 

Using the orthonormal relation (25) it is easy to show that 

/ drrG m (r,r,E) = £ _J__. (31) 
Employing the mathematical formula 

P- - m5(x) (32) 



x + ie x 

and taking the imaginary part of the above equation we have 

Im f drrG m (r,r,E) = —ir ^ 5(E — E mK ). (33) 

Integrating this equation over E from — oo to CT yields 

Im J° dE J drrG m (r,r,E) = -n m n (34) 

where n~ is the number of bound states with negative energies and with angular momen- 
tum mh(when m^Owe have the same number of bound states with angular momentum 
— rah as well). The possibility of a zero-energy bound state will be discussed in the next 
section. In a similar way one can show that 

Im J° dE J drrG$(r,r,E) = 0. (35) 



Here and in (34) the integration over E is performed to the upper limit instead of 
such that it suffers no ambiguity. Combining (34) and (35) we have 

Im f° dE [ drr[G m (r,r,E) - G$}(r,r,E)] = -u^tc. (36) 

On the other hand, substituting (27) and (29) into the rhs of (30) we have 

/ drr[G m (r,r,E)-G^(r,r,E)} = V ^^^^ 77 - (37) 
J na (E - EXk + te) (E - E ma + it) 

Using the radial equations for u$ K (r) and u ma (r) and the boundary condition that these 
radial wave functions vanish at r = and r = R, it is not difficult to show that 

( M mL V u mo) — {E mcr — E(J K )(u(J K ,u ma ). (38) 

Substituting this result into (37) and taking the imaginary part we have 

Im / drr[G m (r,r,E)-G^(r,r,E)} ^ - E^\) -5(E - E ma )]\(u^,u ma )\ 2 . (39) 

Integrating this equation over E from — oo to +oo it is easy to find that 

Im ( +Q ° dE j drr[G m (r,r,E) - G { °\r,r, E)} = 0. (40) 

Similar to the three-dimensional case, this equation means that the total number of states 
in a specific angular momentum channel is not altered by an attractive field, except that 
some scattering states are pulled down into the bound-state region. This result, together 
with (36), leads to 

hnJ*°°dE J drr[G m (r,r,E)-G<®(r,r,E)]=n-7r. (41) 

We have thereupon finished the first step in our establishment of the Levinson theorem. 

The next step is to calculate the lhs of (41) in another way. In the above treatment 
we have discretized the continuous spectrum of E$ K and the continuous part of E mK . In 
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the following we will directly deal with these continuous spectra. We will denote u$ K (r) 
by u^ k (r), and u mK (r) with continuous k by u mk (r), whereas those u mK (r) with discrete 
^(bound states) will be denoted by the original notation. The notations E$ K and E mK 
with continuous k will also be changed to E^ k and E mk ihoth. are equal to h 2 k 2 /2jj) 
respectively. The orthonormal relation for u^ k (r) now takes the form 

(uffi,«20 =*(*-*')■ (42) 

The orthonormal relation for u m k(r) is similar, while that for u mK (r) has the same ap- 
pearance as (25). It is easy to show that 



satisfies the radial equation with V = and the orthonormal relation (42). Thus u^ k {r) 
has the asymptotic form 

(0) / x r-+oo 1 2 I \ 17171 71 



corresponding to (8). In an external field, the wave functions are distorted and thus the 
asymptotic form for u m k(r) becomes 



. ; — x / 2 

u m k{r) — ► W-cos 



7T 



kr- — - - + r] m (k) 



(44) 



corresponding to (9). Note that the coefficient in the asymptotic form is the same for 
u< mk( r ) an d u mk( r )- In this treatment it can be shown that 

G m (r, r>, E)=Y Jw(OO') + [ dk J^X^') (45) 

1 ' ^V^(E-E mK + ie) J V7?(E-E mk + ie) 1 ' 

and thus 

Im / drrG m {r,r,E) = -7T^25{E - E mK ) - it / dk5(E - E mk )(u mk ,u mk ). (46) 
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Integrating this equation over E from to +00 and taking into account the fact that 
EmK < while E mk > 0, we have 



r+00 r i- 

Im dE J drrG m (r,r,E) = -tt J dk (u mk ,u mk ). (47) 

There is no ambiguity in the integration over E since the lower limit is set to be 0~ instead 
of 0(note that E mk may equal 0). As before, the possibility of a zero-energy bound state 
will be discussed in §5. In the same way, we have 

l™[_°°dE J drrG^(r,r,E) = -vr J dk{u<®,u®). (48) 

It should be pointed out that the integrands in both (47) and (48) are 5(0) (= 00) according 
to Eq.(42) and a similar equation for u mk (r). However, there is a subtle difference between 
these two 5 functions, and it is this difference that leads to the Levinson theorem. Since 
both integrands are singular, we first evaluate 

(u mk ,u ml ) ro - (uS> u 2!)ro = / dru *mk( r )u m i(r) - [ dv u®£ (r)u^j(r) ( 49 ) 

Jo Jo 

where r is a large but finite radius, and finally take the limit / — > k and r — > 00. As in 
the discretized case, we have the boundary condition 

«S(0) = 0, Umk(0)=0. (50) 
Using the radial equation and this boundary condition it is easy to show that 

{k 2 - l 2 ){u mk ,u m i) ro = u* mk {r )u' ml {r ) - u* k (r )u ml (r ) . (51) 
Since r is large, we can use (44) to evaluate the rhs and in the limit I — > k we get 

(umk,u mk ) rQ = — + -r]' m (k) - — — cos[2/cr + 2r} m (k)]. (52) 

71 71 Z7TK 
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In the same way we have 



(«2,«a) rt> = ^-^oo8 2AT . (53) 



Therefore, 



(umk,u m k)ro ~ («S> u ml)r = -v'm( k ) H 7r~ sm2i] m (k) H cos 2krQ sin. 2 1] m (k) 

71 Z TTK 

(54) 

where we have employed the well-known formula 

hm = 5(h). 

So far in this section V(r) need not be V a (r). In the following we set V(r) = V a (r). Then 
Eq.(18) is valid, and we have sin2?7 m (0) = 0. Therefore the second term on the rhs of 
(54) vanishes. Integrating this result over k (from to +oo), taking the limit r — > oo, 
and incorporating the results (47), (48) we arrive at 

lm£°° dE J drr[G m (r,r,E) -G$(r,r,E)] 

= Vm(0)-Vm(oc) - jim jT^dfc^^siiiVW- (55) 

The last term in this equation can be decomposed into two integrals, the first from to 
e = + , which can be shown to vanish on account of (18), while the second from e to +oo, 
which also vanishes in the limit r — > oo since the factor cos 2kr oscillates very rapidly. 
Therefore we have 

hn J*°°dE J drr[G m (r,r,E)-G%\r,r,E)]=r tm (p)-r tm (oo). (56) 
Combining (56) and (41) we arrive at the Levinson theorem: 

Vm(0) - ^m(oo) = V, m = 0,1,2,.... (57) 
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In the next section we will discuss some aspects of this theorem, and get a more general 
form which takes zero-energy bound states into consideration. 

§5. Discussions 

In this section we discuss and clarify some aspects of the Levinson theorem obtained in 
the form (57) in the last section. 

1. On zero-energy bound states. In §4 we have not taken into account the possible 
existence of a zero-energy bound state. Indeed, this may occur for a square well potential 
with radius a and depth Vq when m > 1 and J m _i(£) = where £ = k^a and ho = 
y/2fjVo/h. (For m — 0, 1 regular solutions with zero energy may be found but they are 
not normalizable and thus are not bound states, see below) The existence of a zero-energy 
bound state would not alter the results (34)- (36) where n~ is the number of bound states 
with negative energies. Therefore Eq.(41) remains valid in this case. On the other hand, 
Eq.(47) becomes 

r+oo r r 

Im dE J drrG m (r,r,E) = -n - n J dk{u mk ,u rnk ) (58) 
which is an obvious consequence of (46). Accordingly, the result (56) becomes 

Im J^dE J drr[G m (r,r,E)-G<V(r,r,E)] = Vm (0) - Vm (oo) - n. (59) 
Hence, the Levinson theorem takes in the present case the form 

Vm(0) - Vm(00) = n m 7T (60) 

where n m = n~ + 1 represents the total number of bound states, including the one with 
zero energy. When there is no zero-energy bound state, n m = n~. Thus Eq.(60) holds in 
any case and is the final form of the Levinson theorem. 
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To see the difference in zero-energy states between two and three dimensions, we notice 
that the two-dimensional radial wave function u mK (r) satisfies (24). When regarded as a 
one-dimensional Schrodinger equation, the effective potential reads 

f „, ,„ , A 2 (m 2 -l/4) 

where the subscript "2" indicates two dimensions. In three dimensions, with ipi mK (r, 9, <fi) = 
r ~ 1 XiKij')Yim{0 , 0)(here r, 6 etc. should not be confused with those in two dimensions), 
the radial wave function Xik(j') satisfies 



x'L + 



h 2 



XiK = 0, / = 0,1,2,.... (24') 



is 

V 3 (r) =V(r) + 



When this is regarded as a one-dimensional Schrodinger equation, the effective potential 

% 2 l(l + 1) 
2/xr 2 

which is obviously different from V 2 (r), and as a consequence the zero-energy solutions 
are different from those of (24). More specifically, with V(r) = V a (r) and E m0 = 0, the 
exterior solution(r > a) of (24) reads 

<o(r) = r- m+1 / 2 . 

A zero-energy solution exists if V a (r) is such that the interior solution u m0 (r) can be 
connected with «^ (r) at r = a smoothly. This leads to the above mentioned condition for 
a square well potential. Here we are concerned about the normalizability of the solution. 
It is clear that the above solution can be normalized and thus is a bound state only when 
m > 1. We have ipi r —> 0, and ^oo is finite at infinity though u 00 (r) oo, thus both 
are regular. But either of them can not be normalized and thus is not a bound state. 
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In this case r] m (0) gets no additional term and Eq.(57) is not modified. When m > 1, 
however, n~ — > n~ + 1 = n m and r) m (0) gets an additional n if a zero-energy solution(a 
bound state) actually exists. In three dimensions, with V(r) = V a (r) and E m = 0, the 
exterior solution of (24') reads 

xU r ) = r ~ l - 

This is normalizable when / > 0. Thus only / = is distinguished. When the / = solution 
really emerges, 5o(0) gets an additional 7r/2 as pointed out in §3, and the Levinson theorem 
is modified by n — > n + 1/2. This is different from the case m — 0, 1 in two dimensions 
where ?7 m (0) always equals a multiple of ir. When / > the case is basically the same as 
m > 1 in two dimensions. 

2. About r) m (oo). Write down two radial equations for u mk (r) and u mk {r) m the 
external potentials U(r) and U(r) respectively. Using the boundary condition (50) and 
the asymptotic form (44) for u m k(r) and similar ones for u m kij'), h is easy to show that 

sin[?7 m (£;) - fj m (k)} = / dr [U (r) - U (r)]u* mk (r)u mk (r) . (61) 

n k Jo 

Now we take U(r) = U(r)(not necessarily be V a (r)) and U(r) = 0, it is natural to define 
Vm(k) = in the absence of an external field. Obviously, u m k{ r ) = u mk( r )> so (61) becomes 

sin r] m (k) = f dru* mk (r)V(r)u^ k (r). (62) 

For very large k, V(r) can be ignored in the radial equation since it is less singular 
than r~ 2 near the origin and well behaved elsewhere as assumed and so can be neglected 
everywhere (cf. Eq.(24)). Therefore u mk (r) in (62) can be replaced by u^ k (r) in this limit, 
and we have for very large k 

sin Vm (k) = -^fdr\u { °l(r)\ 2 V(r). (63) 
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Substituting u^ k (r) = \^krJ m (kr) into this equation, and using the asymptotic formula 
for the Bessel functions at large argument, we have approximately 



Further, since k is very large, the cosine oscillates very rapidly and thus the squared cosine 
may be replaced by its mean value 1/2. In this way we arrive at 



for very large k, where we assume that the integral exists. This result has the same 
form as that in three dimensions. Of course it holds in the special case V(r) = V a (r). 
Obviously, sinr] m (k) — > when k — > oo, hence we can freely define ?7 m (oo) = 0. Under 
this convention the Levinson theorem takes the form 



This means that the phase shift at threshold serves as a counter for the bound states. This 
is similar to the case in three dimensions, but is somewhat simpler. In three dimensions, 
the case / = should be modified when there exists a zero-energy resonance(a half-bound 
state). Here we have no such problem. 

3. Extension to more general potentials. In the above development of the Levinson 
theorem, we have assumed that V(r) = when r > a. We also assumed that V(r) is less 
singular than r~ 2 when r — > 0, so that the power dependence in (13) near r = is valid. 
Indeed, the existence of the integral in (65) requires that V(r) is less singular than r -1 
when r — > 0. We will always assume that V(r) is sufficiently well behaved near the origin 
such that the above requirements are all satisfied. On the other hand, the radius a beyond 
which V(r) vanishes is not specified in our discussion. Though both r] m (0) and n m in (66) 




(64) 




(65) 



^ m (0)=n m vr, m = 0,1,2,.... 



(66) 
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depend on the particular form of V(r) and thus depend on a, the equality between them 
does not. Hence one expects that (66) remains valid when V a (r) is varied continuously to 
the limit a — > oo if n m remains finite in the process. It seems that the Levinson theorem 
holds for quite general potentials(well behaved near the origin as emphasized above) as 
long as they decrease rapidly enough when r — > oo such that the total number of bound 
states in a particular angular momentum channel is finite. 

Extension of the present work to relativistic quantum mechanics is currently under 
progress. 
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